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Abstract 



We calculate analytically the highly damped quasinormal mode spectra of generic single- 
horizon black holes using the rigorous WKB techniques of Andersson and Howls ^B]- We 
thereby provide a firm foundation for previous analysis, and point out some of their possible 
limitations. The numerical coefficient in the real part of the highly damped frequency is 
generically determined by the behavior of coupling of the perturbation to the gravitational 
field near the origin, as expressed in tortoise coordinates. This fact makes it difficult to 
understand how the famous ln(3) could be related to the quantum gravitational microstates 
near the horizon. 



1. Introduction 

The highly damped limit of black hole quasinormal modes (QNMs) has recently at- 
tracted a great deal of attention, in large part due to the work of Hod£Q and Dreyer[2]. 
These authors used different, but closely related, arguments to extract information about 
the underlying black hole quantum gravitational microstates from the real part of the highly 
damped QNMs. Their arguments can be summarized loosely as follows: For Schwarzschild 
black holes the real QNM frequency seems to approach a unique value (independent of 
the nature and angular momentum of the perturbation) as the imaginary part goes to 
infinity [3]. From this fact, one might conclude, as did Hod 1., that this limiting value 
constitutes a fundamental oscillation frequency associated with the dynamics of the event 
horizon. Dimensional arguments require this frequency to be proportional to the Hawking 
temperature: 

k 

to = a ^T BH , (1.1) 

where a is a dimensionless constant 1 . 

Given the existence of such a frequency, semi-classical arguments require the existence 
of states in the energy spectrum that are separated by the corresponding energy quantum 

AE n = huoAn , (1.2) 

where n is the integer labeling the states and An = 1. In the large n limit this expression 
can be integrated to yield 

rdE 1 r dE 1 
n = — = - / 7jr- = -Sbh ■ (1.3) 
J lu a J 1 BH Oi 

By the first law of black hole thermodynamics, this is proportional to the Bekenstein- 
HawkingjlJ |5| entropy, which therefore has states that are equally spaced in the semi- 
classical limit: 

Sbh = an . (1.4) 

That the black hole area/entropy is an adiabatic invariant with equally spaced spectrum 
was first conjectured by Bekenstein[6j and later by Mukhanovj7]. The above form of the ar- 
gument applies to all black holes described by a single dimensionful parameter and assumes 
only the existence of a fundamental vibration frequency associated with the horizon. 



■""Note that since the Hawking temperature is proportional to h the right hand side is a purely classical 
quantity, as required. We will henceforth use units in which h and Boltzman's constant k are unity. 
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It has often been claimed that an equally spaced spectrum appears to contradict Hawk- 
ing's prediction that black holes radiate energy with a thermal spectrum. That is, the 
spacing between the resulting emission frequencies are of the same order as the Hawking 
temperature, even for macroscopic black holes to which the Hawking's semi-classical ar- 
guments should in principle apply. In this regard it is important to remember that the 
arguments leading to ()1.4|) imply only that the area spectrum of spherically symmetric 
black holes should contain states that are equally spaced, with separation a. This does not 
necessarily imply that the entire black hole area spectrum is equally spaced. Moreover, 
this prediction applies only to the case of exact spherical symmetry. It is possible, and 
in fact very likely, that the addition of further degrees of freedom (e.g. deviations from 
spherical symmetry), will lead to a sufficient number of new lines in the spectrum so as to 
make it effectively continuous, at least for macroscopic black holes in which the principle 
quantum number n is very large, and the Hawking temperature is very small. 

Hod was the first to notice from numerical calculations [3] that for real part of the 
highly damped QNMs of 4-d Schwarzschild black holes, the constant of proportionality a 
appeared to be ln(3). This leads to a semi-classical Bekenstein-Hawking entropy spectrum 
of ln(3 n ) which has a natural statistical mechanical interpretation in terms of a black hole 
horizon made out of n elements of area, each with three internal degrees of freedom. 

Dreyer j2] discovered that Hod's arguments appeared to be consistent with the statistical 
mechanical black hole entropy that arises in loop quantum gravity [Sj, and remarkably 
could be used to independently and consistently fix the elusive Barbero-Immirzi parameter. 
Unfortunately a recent recalculation [§] of the statistical mechanical entropy of black holes 
in LQG seems to bring into question the validity of Dreyer's analysis. 

These developments subsequently lead to numerous analytic analysis of the highly 
damped QNMs of generic black holes. The first was due to Motl^U] who proved that 
the coefficient a for 4-D Schwarzschild black holes was indeed precisely log(3), verifying 
the conjecture by HodpQ. Motl and Nietzke^J were able to prove that the ln(3) per- 
sists for Schwarzschild in all dimensions above three, as argued in [T2]. Birmingham and 
Carlip[13] then showed the existence of an intriguing, but somewhat different, relationship 
between the highly damped QNMs of the BTZ black hole and its quantum gravitational 
states. There are by now many analytic calculations of the highly damped QNM spectrum 
for a large variety of black holes|14|. It seems that for spacetime dimensions 4 and up, the 
intriguing conjectures described above apply at most to single horizon, asymptotically flat 
black holes. 
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Most calculations of the highly damped QNMs focus on particular black hole solutions. 
This includes a recent paper by Tamaki et a/J3], which start with a general metric,but 
then specialized to 4-D Schwarzschild and dilaton black holes. More recently, two papers 
have calculated the highly damped QNMs for large classes of single horizon black holes 
within a single formalism. Kettner et al[16\ analyzed non-minimally coupled scalar fields 
in the background of black holes in generic 2-d dilaton gravity. In an elegant analysis, 
Das et looked at the QNMs of scalar fields in the background of completely general 

static, spherically symmetric black holes in d dimensions. It turns out that the two classes 
considered in these papers are closely related. 

The purpose of the present paper is to re-examine the most general class of spherically 
symmetric, asymptotically flat black holes using the rigorous WKB techniques of Andersson 
and Howls ^H]. We are thereby able to provide a firm foundation for previous analysis, 
and point out some of their possible limitations. Our general formalism also sheds light 
on the source of the universality and proposed physical significance of the famous ln(3) 
in the real part of the highly damped QNM frequency. We show in fact that the ln(3) is 
not generic in principle, despite the fact that in practice it does appear valid for higher 
dimensional, asymptotically flat black holes even in the presence of other dynamical fields 
(i.e. not just for higher dimensional Schwarzschild). Moreover, the value of this coefficient is 
generically determined by the exponent that determines the behavior of the coupling of the 
perturbation to the gravitational field near the origin, as expressed in tortoise coordinates. 
Although there may be general arguments that restrict the value of this parameter, it is 
hard to see how it could be related to the quantum gravitational properties of the horizon. 

In Sec. 2, we set up the problem by introducing a general equation for scalar pertur- 
bations in a static, spherically symmetric black hole background. In Section 3 we present 
some general properties of the large damping limit of this equation, and make some general 
observations about the universality of the resulting frequencies. In Sec. 4, we investigate 
the global structure of the Stokes and anti-Stokes lines for the highly damped QNMs of a 
specific, but still very large, class of black holes spacetimes. In Sec. 5, we calculate the 
corresponding QNM frequencies using the WKB techniques presented in [TH]. In Sec. 6, 
we investigate other possibilities for the change of variable, which is introduced in Sec. 4 
to obtain suitable behavior of the Stokes and anti-Stokes lines. A summary and discussion 
is given in Sec. 7. 
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2. General Formalism 

The general form of the wave equation that we wish to consider can be suggestively 
written in 2-dimensional form 



where g^ v is the two metric on the (r, t) plane, and is a scalar with respect to 2-d 
coordinate transformations. Both the metric and dilaton are assumed static, so that one 
can find coordinates (x, t) in which <p = <j>(x) and the metric takes the form 

ds 2 = -f{x)dt 2 + -^—dx 2 . (2.2) 

Note that we are as yet making no assumptions about the dynamics that give rise to 
this metric. The functions f(x), g(x), and h(x) = h[<p(x)] are completely arbitrary. By 
further restricting the coordinate system, it is possible to eliminate at most one of these 
functions, so the system is in fact completely specified by two arbitrary functions. V(x) 
is determined by the properties of the scalar field. For higher dimensional black holes it 
contains the centripetal term associated with the angular momentum of the perturbation. 

In order to restrict to single horizon black hole spacetimes we assume that h(x) is 
monotonic and vanishes at x — 0, which is a singular point in the spacetime. Moreover, 
we assume f(x) and g(x) have simple zeros at the same non-zero Xh, the horizon location. 
Their ratio H(x) = is assumed to be a regular, nowhere vanishing, analytic function 

of xjni. 

At this stage it is useful to make contact with two previous analysis. The first is the 
calculation of Das et al\£j\ 2 in which the most general static, spherically symmetric metric 
in d spacetime dimensions was considered: 

ds 2 = -f(r)dt 2 + -^--dr 2 + r 2 d^ n) , (2.3) 
9{r) 

where n = d — 2 and dil^ is the line element on the unit n-sphere. By dimensionally 
reducing the wave equation for a minimally coupled scalar field in this background and 
making the identifications x = r and h(x) = r n , one obtains precisely (|2.1j) . with 

V = r"'^;- 1 ) (2.4) 



2 The authors are grateful to S. Das for a very useful discussion explaining the significance of these 
results. 
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Alternatively, one can consider (|2.1J1 . as in [TBj . to be the equation describing a scalar field 
non-minimally coupled to a black hole metric and dilaton in generic vacuum 2-d dilaton 
gravity. In this case one can choose <fi = x to find that[P9*] 

f(x)=g{x) = J(x)-2GM, (2.5) 

where M is the mass of the black hole and J(x) is determined by the dilaton potential, 
which is different for different theories. 

These two formalisms coincide in the special case that the metric in (|2.3J1 describes the 
higher dimensional Schwarzschild spacetime, for which 

J(0) oc 0^ 1/n . (2.6) 



An important observation about ()2.1j) is that the 2-dimensional D'Alembertion is in- 
variant with respect to conformal transformations of the metric. In particular, if one writes 

2" 



ds 2 = f{x) 



_ o ( dx 
-dt 2 + 



(2.7) 



F(x) J 

the left hand side of the wave equation depends only on the coupling h(x) and F(x) = 
f(x)g(x): 



- ^r\d 2 ^(x } t) + d x [h(x)F(x)d x iP(x, t)] = J f -V{x)^{x, t) . (2.8) 
F(x) V 9 

In terms of the redefined field 

^ t )=*M == e*> \ (2.9) 
yjF(x)h(x) 

one obtains 

d 2 ^ 

— + ii!(*)* = 0, (2.10) 

where 

R i x ) = ^T-^) > (2- 11 ) 
t z {x) 

and the scattering potential between the dilaton and the matter perturbation is given by 



/N IF" 1 ( F'V Ih" lF'h! ifh'V J{V{x) 
with prime denoting differentiation with respect to x. 
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Equation (j2.1(jj) has singular points at the horizon Xh, where F(x) vanishes but h(x) is 
assumed regular, and at the origin x — 0, where h(x) vanishes. Furthermore, the potential 
U(x) is assumed to vanish as x — > oo faster than 1/F 2 . 

For future reference we note that in terms of the tortoise coordinate 

AT 

F[x) 



and a rescaled wave function \1/ = / y/F, the wave equation (|2.1U|) takes the simple form 

dz- 

where 



+ (aj 2 -U h (z))y = , (2.14) 



1 h" 1 (h'V F , 

and the prime here denotes differentiation with respect to z. Under the assumptions we 
made for F(x), h(x), and U(x), the potential Uh goes to zero at both the horizon (z — > — oo) 
and spatial infinity (z — > oo). Therefore the asymptotic behavior of the solutions in those 
regions are 

try/ - \ f e^ z asz^ -oo (x -> ar fc ) , , , 

*W « i M . oo (x - oo) . (2 ' 16) 



3. Large Damping Limit and Universality 

Equation ()2.10j) with boundary conditions ()2.16|) defines the general problem of finding 
QNMs for a generic black hole. We now consider the QNMs in the large damping limit 

\uj 2 \ -> [Im(cj)] 2 -> oo . (3.1) 

Since we will be using complex analytic techniques, the behavior of U(x) on the entire 
complex plane may be relevant. However, in the large damping limit case, the uj 2 /F 2 term 
in R(x) will dominate U(x) everywhere, unless one of the terms in U(x) diverges. According 
to our assumptions, this can only happen at the origin, or at the horizon. However, since 
the uo/F 2 also diverges at the horizon, it will dominate there as well in the large damping 
limit. Thus, U(x) is essentially irrelevant in this limit except near the origin. Suppose 
without loss of generality that h — > x a , F — > ar, and V — > x 1 in this limit. As long as 
a + (3 — 7 < 2, which is the case for all previously considered calculations, V(x) can be 
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(3.3) 



neglected. Thus, for the highly damped QNMs, R(x) can be approximated on the entire 
complex plane by 

RM ^ (<» + /?)[(<» + /?) -2] 

Note that in the WKB analysis where the two solutions to Eq. (j2.1()j) are 

' ^f{ x ) = Q- l / 2 ( x ) exp [+i J t x Q{x')dx'\ , 

= Q- l ' 2 {x) exp [-i J t x Q{x')dx'\ , 

we need to deal with Q rather than R. Here Q 2 = R+ extra term, where the extra term is 
to make sure that the WKB solutions (J3.3)) have the appropriate behavior near the origin 
of the complex x-plane. From Eq. (|3.2j) we have R ~ —(a + (3)[(a + 0) — 2]/4x 2 when 
x — > 0. Using Eq. (j2.1()|) . it is easy to show that \I> ~ a; i/2±(a+/3-i)/2 ag x _^ q Therefore, 
we need to take Q 2 = R — l/4x 2 , so that the WKB solutions ^2 have the correct behavior 
of the form \1/ ~ a; 1 / 2 ±( a +/ 3 - 1 )/ 2 close to the origin. As a result, the approximate behavior 
of Q 2 on the entire complex plane is 

=^)"^~^ ^^2^ • ( 3 - 4 ) 

Equation (|3.4j) is the key to understanding the universality of the highly damped QNM 
spectrum for generic single horizon black holes. The analytic WKB techniques used in the 
present paper allow the QNMs to be determined purely by the structure of the poles at 
x = and x = Xh- In fact, as we will show, when the calculation can be done using these 
techniques then the QNMs are determined by the relationship 

e^r = -(l + 2cos((o-l)7r)) , (3.5) 

where 



IF 

is the surface gravity of the black hole and have defined the parameter 



a 



(3.6) 



(3.7) 



(1-/5) 

The parameter a is in fact the precise combination of a and f3 that is invariant under 
coordinate transformations of the form x — >• x = x q . In fact, it is easy to verify that a gives 
the rate at which the coupling h(x) approaches zero as a function of the tortoise coordinate 
z. That is: 

h[x{z)} -> z a (3.8) 
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as x — > 0. 

At a heuristic level, one can see that only a and the surface gravity can affect the values 
of the highly damped QNMs by noting that in tortoise coordinates Uh — > ■ Thus for 
large |u;| 2 , the equation has a universal form with different black hole backgrounds being 
distinguished only by the residues of the poles at x = and x = Xh, i.e. a and k respectively. 
Moreover, the parameter a is determined by the coupling of the perturbation to the metric, 
and does not seem to have any direct relationship to properties of the horizon. 

We close this section by noting that a = 1 for scalar fields in the background of the 
Schwarzschild metric in D = n + 2 spacetime dimensions. Since a is coordinate invariant, 
we will work with the radial coordinate r. In this case: 

F(r) = (1 - 2GM/r D - 3 ) (3.9) 

where M is the ADM mass as measured at infinity. Moreover, h(r) = r D ~ 2 so that 
a = D-2tmd(3 = 3-Dto give a = (D - 2)/(l - (3 - D)) = 1. This in turn leads to 
the known result 

^— = (2n + 1)ttz + Zn(3) . (3.10) 
Tbh 

It is interesting to note that a = 1 also for stringy black holes in 4 and 5 dimensions (see 
|17j . in which a = qD/2). It seems important to understand the necessary and sufficient 
conditions on asymptotically flat single horizon black holes which lead to a = 1. These 
conditions, coupled with the analysis given above, would account for the universality of 
the highly damped corresponding QNMs for single horizon black holes. They would also 
presumably give further insight as to whether the ln(3) can be related to the underlying 
quantum gravity theory. 

4. 2-d dilaton gravity with power-law potentials 

So far we have kept the arguments as general as possible. For concreteness, we now 
focus on a generic class of black holes in 2-dimensional dilatonic gravity. In the conclusions 
we explain the circumstances under which the techniques and results in the present section 
can be applied to a much more general setting. 

With the gauge choice of (j) — x, and the "power-law potentials" considered in [16,, we 
have 

J(x)=x 1 ~ b . (4.1) 
9 



The power-law metric is 

ds 2 = -f{x)dt 2 + ^f- , (4.2) 

where 

fix) = J(x) - 2GM , (4.3) 
which locates the black hole horizon at 

x h = [2GM) 1/1 - b . (4.4) 

We further need to restrict consideration to the sub-class of "power-law potentials" for 
which 

< b < 1 . (4.5) 

The upper bound is necessary for the existence of black hole solutions, while the lower 
bound guarantees that the physically relevant solutions of the action of the 2-d dilatonic 
gravity are asymptotically flat. 

We also need to know the behavior of the coupling parameter h(x) between the dilaton 
and matter field in the complex x-plane. In the physically motivated case of spherically 
reduced Einstein gravity the natural choice for the coupling h is h = xj201- Therefore it 
would be reasonable to use h = x a , where a is an arbitrary parameter. Substituting h = x a 
and V(x) = in the scattering potential (|2.12j) . and keeping only the terms which diverge 
faster — » we get 

U(x) = \^ . (4.6) 

Note that (3 = in this parameterization. 

To calculate the QNM frequencies in the highly damped limit where Im(u;) — > — oo, we 
will follow the analytical method based on WKB approximation suggested by Andersson 
and Howls in Ref . JH] , which ultimately provides the motivation for the somewhat simpler, 
monodromy method of Refs. [10] and [TT] that was applied in Refs. (THj and [T7j . 

The zeros and poles of the function Q in the WKB solutions ()3.3|) play a crucial role 
in determining the behavior of the Stokes lines on which / Qdx is purely imaginary and 
anti-Stokes lines on which / Qdx is purely real. These lines are vital for the WKB analysis 
of Eq. (|2.10j) . From each zero three Stokes lines and three anti-Stokes lines emanate. The 
angle between Stokes (anti-Stokes) lines is 120°, and the angle between Stokes and anti- 
Stokes lines is 60°. The poles of the function Q 2 are located at the origin x = 0, and the 
event horizon x = Xh- It is not difficult to show that the zeros of Q 2 approach the origin of 



10 



Figure 1: A schematic illustration of Stokes (dashed) and anti-Stokes (solid) lines in the 
complex x-plane. The hollow circles represent the two zeros of Q 2 , while the filled circles 
are the two poles at the origin and the event horizon. 



the complex x-plane when Im(a>) — ► — oo. In this limit, Q 2 may be simplified by expanding 
Eq. (j2.11j) . with U(x) to be the scattering potential ()4.6)1 . in a power series near x = 0; 

q* = r _JL K ^!_(^i£. (4.7) 

Note that since in this model h = x a , and f(x) approaches to a constant when x — ► 
we can derive this equation by simply substituting /3 = in Eq. (|3.4jl . It is clear from 
Eq. ()4.7|) that the function Q 2 has only two simple zeros located on the imaginary axis of 
the complex x-plane. The schematic behavior of the Stokes and anti-Stokes lines on the 
entire complex plane is displayed in Fig. ^ As one can see, there are no anti-Stokes lines 
extending to infinity. The unbounded anti-Stokes lines are crucial ingredients in solving 
for QNM frequencies in the analytical methods suggested in Refs. ^U], ^1] and ^H]- We 
have found no way of solving this problem without the presence of such anti-Stokes lines. 
In order to overcome this problem we use a change of variable of the form 

x = y q . (4.8) 

In this new coordinate, the tortoise-like (spatial) coordinate is defined by 

dz = ^ , (4.9) 

f{y) 



n 




Figure 2: A schematic illustration of Stokes (dashed) and anti-Stokes (solid) lines for 
b = 1/5, 2/7, and 3/4 with q = 5, 7, and 4 respectively in the complex ?/-plane. The hollow 
circles represent the zeros of Q 2 , while the filled circles are the poles at the origin, the 
event horizon, and the fictitious horizons. For to = lor + iuj, with ujr > and ujj < 0, the 
anti-Stokes line labeled e asymptotes to a straight line at an angle arctan(— ui /ujr) / qb. 



m = — • (4-io) 



where 

2,9(1-6) _ 2GM 
qyi 

The new horizon is located at 

y h = {2GM) 1/q[1 - b) . (4.11) 

Since f(y) — > y 1 ~ q /q as y — > 0, and ft = y 9 " in this new coordinate, the behavior of the 
function Q 2 close to the origin of the complex y-plane can be derived by simply substituting 
1 — (3 — q, and a = qa in Eq. (|3.4|) . where we get 



2 -R 1 q2{a ~ l)2 (412) 

Q - R -^2~p- V " ( 4 - 12 ) 



We now have the freedom to increase the number of zeros of the function Q 2 with the help 
of our arbitrary parameter q. The important point is to choose q in such a way that we have 
enough number of zeros to accommodate not only for the anti-Stokes line encircling the 
event horizon on the positive real axis, but also for anti-Stokes lines extending to infinity 
in the complex y-plane. For a particular value of b in the region < b < 1, there may be 
more than one choice of the parameter q that produces the desired behavior of anti-Stokes 
lines. To make matters concise however, let us assume that 6 is a rational number which 
may be written as the ratio of two integers 

b = - , (4.13) 

n 
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where m < n. We then conjecture that the choice of q = n will always produce the desired 
structure of anti-Stokes lines for any rational b in the interval < b < 1. The behavior 
of Stokes and anti-Stokes lines are schematically plotted in Fig. [21 for particular values of 
b = 1/5, 2/7, and 3/4 with q = 5, 7, and 4 respectively. Note that in addition to the event 
horizon on the positive real axis, we may also have extra fictitious horizons in the complex 
y-plane. The horizons are located at 

where k = 0, 1, 2, n — m — 1 with k = corresponding to the event horizon. Since the 
number of zeros of the function Q 2 is 2q = 2n, the anti-Stokes lines labeled e and d in 
Fig. El emanate with an angle 3ir/n with respect to each other. In addition, according to 
Eq. (J4.14j) . the argument of the first fictitious horizon in the complex y-pl&ne is 2ir/(n — m). 
This means that the angle between the two fictitious horizons on either side of the event 
horizon is 47r/(n — m). Based on our considerations 3ir/n < 47r/(n — m). This means that 
our fictitious horizons will always stay to the left of the anti-Stokes lines labeled d and e 
in Fig. H 

The behavior of the Stokes lines on which Re(/ Qdy) = and anti-Stokes lines on which 
Im(J Qdy) = near the origin of the complex y-plane may be determined using Eq. (|4.12|) . 
Let us investigate the behavior of these lines away from the origin. In this region we may 
write 



qy q 1 lu \ 2 ( ny n 1 



to 



Q ~ \y^-b) - 2GM) [yn-™-2GMJ ' (4 ' 15) 
when we take b = m/n and q = n. To determine the behavior of the Stokes and anti-Stokes 
lines we need to evaluate / Qdy. Integrating Q with respect to y gives 



n „ n n ~ m ~ 1 



Qdy~±{-y™ + — £ y k 



fc=0 



m-2 m-l-p / „. n 

E (y^r—, — + ivkT- 1 Mi- 1 

t^o m-l-p V Vk, 



(4.16) 

where we added a constant to make sure that / Qdy = at y = 0. It is clear that / Qdy is 
a multivalued function. To make this function single valued we need to introduce branch 
cuts in the complex y-plane emanating from each of the horizons. Near the horizons this 
equation may be simplified to give 



/ 



Qdy~± 



-^rm [1-3- 

n-m V y k 



u . (4.17) 
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Assuming u to be purely imaginary, it is easy to see that Im(J Qdy) is single valued around 
Uk with k = and k — (n — m)/2 (if n — m is even), but around all the other fictitious 
horizons, where y^ is complex, Im(J Qdy) is multivalued. This means that anti-Stokes lines, 
on which Im(J Qdy) = 0, are well defined around the horizons located on the real axis and 
therefore are allowed to cross the branch cuts emanating from these horizons. Around all 
the other fictitious horizons the anti-Stokes lines are not well defined and therefore are not 
allowed to cross the branch cuts. Some of the anti-Stokes lines may terminate on these 
fictitious horizons or on the branch cuts emanating from them. It is also easy to see that 
Re(J Qdy) is multivalued around all the horizons except around y^ with k = (n — m)/4 and 
k = 3(n — m)/4 (if n — m is even). Therefore Stokes lines, on which Re(J Qdy) = 0, are 
not allowed to cross the branch cuts emanating from the horizons which are not located 
on the imaginary axis. 

As we mentioned earlier, the fictitious horizons are always located to the left of the 
unbounded anti-Stokes lines d and e. Therefore we are able to introduce the branch cuts 
emanating from the fictitious horizons in a way that they do not affect the behavior of the 
Stokes and anti-Stokes lines on the right hand side of the complex y-plane. We will see 
in the next section that the behavior of the Stokes and anti-Stokes lines to the left of the 
lines d and e are irrelevant in our calculations and therefore have not been plotted in Fig. 
121 For the interested reader we plot the behavior of the Stokes and anti-Stokes lines on the 
entire complex y-plane for 6 = 3/7 with q = 7 in Fig. El Note that in this figure two of the 
anti-Stokes lines terminate on the fictitious horizons located on the imaginary axis. For 
more examples on the behavior of anti-Stokes lines in the vicinity of the fictitious horizons 
refer to Ref . [H] . 

5. Quasinormal mode frequencies 

In the previous section we conjectured that with the choice of q = n we may always 
have a desired structure of Stokes and anti-Stokes lines similar to the diagrams shown in 
Fig. |21 for any rational b = m/n . In this section we will show how to solve for the highly 
damped QNM frequencies using the diagrams in Fig. [21 

We start at point e on the anti-Stokes line extending to infinity with an outgoing wave 
solution 

^ e = ^S tl} . (5.1) 
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Figure 3: A schematic illustration of Stokes (dashed) and anti-Stokes (solid) lines for 
6 = 3/7 with q = 7 in the complex y-plane. The hollow circles represent the zeros of 
Q 2 , while the filled circles are the poles at the origin and the horizons. The branch cuts 
emanating from the horizons are illustrated with zigzagged lines. 

Note that in order to have the outgoing-wave solution at infinity to be proportional to 
and the ingoing- wave solution at the event horizon to be proportional to ^2, we have 
chosen the phase of the square-root of Q 2 such that Q = \J R— l/4y 2 ~ qy qh ~ X LO as y — > 00. 
To do this, we need to introduce branch cuts emanating from the simple zeros of Q 2 . These 
branch cuts may usually be placed in such a way that they do not affect our analysis. 

To go from e to / we cross a Stokes line emanating from t\ on which ^ is dominant 
according to our choice of phase for Q. When we cross a Stokes line we need to account 
for "Stokes phenomenon" [21] which means that the coefficient of the dominant solution 
(\l/i in this case) on the relevant Stokes line remains unchanged and the coefficient of 
the subdominant solution i n this case) picks up a contribution proportional to the 
coefficient of the dominant solution. The constant of proportionality is called "Stokes 
constant" . In our case where we are dealing with isolated simple zeros, the Stokes constant 
is simply +i for crossing in the positive (anti-clockwise) direction, and —i for crossing in the 
negative (clockwise) direction. Therefore, after crossing the first Stokes line, our solution 
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becomes 



(5.2) 



This solution will not change in character when we move along the anti-Stokes line labeled 
/ from t\ to £2 because no Stokes lines cross this contour, but we need to change the lower 
limit of the phase-integral of ^1 and ^2 to £2- This means that we need to evaluate the 
integrals of the type 



Qdy 



^fo- 1 ) 2 q 2 {a-l) 2 



q'y 



11/2 



-LU 



dy 



(5.3) 



(2GM) 2 Ay 2 

where £j and tj are two simple zeros of the function Q 2 . We can solve this integral by 
introducing a new variable 77 = 2ujy q / [2G M (a — 1)] which maps the zeros to —1 or +1, and 
this integral becomes 

\ 1/2 

>2 



a(a 



+ 1 f ±i I 1 \*/- a _i 
7« = J o / . ( 1 " - ) *7 = T^vr . (5.4) 



2 J^i \ rj 2 J 2 

We now can change the lower limit of the phase integral in ()5.2j) to £2 by simply writing 



(*2) 



-«7i2 



(*2) 



(5.5) 



We may extend this solution to g by crossing another Stokes line on which is dominant. 
Hence, at g we obtain 



y g = e ni2^(*2) _ ■ (y 7l2 + e -i 712 -j ^ 



(to) 

2 



(5.6) 



We then proceed from (7 to g' on the anti-Stokes line that loops around the pole at y^. 
This contour does not cross any Stokes line, but we need to change the lower limit of the 
phase integral to t 3 . To do this, we may write 



o«7l2^723 



(ts) 



i (V 712 + e" i712 ) e -^ 23 ^ ( * 3) 



2 5 



(5.7) 



where 723 = T + 723. Here T is the integral of Q along a contour encircling the pole at y^ 
in the clockwise direction and may be evaluated using the fact that 

1 



r = f Qdy = —2ixi Res Q = —2ni 
J " y=\yh\ 



-(2GM) 



6/(1-6) 



In terms of the generic black hole Hawking temperature 



T, 



K 



BH 



2tt Att dx 
16 



(5.8) 



(5.9) 



r is therefore 



T = -i- 



10 



2T, 



BH 



(5.10) 



We now connect the solution to the point h by stepping inside the anti-Stokes line 
encircling the horizon. We cross the first Stokes line on which \l/ 2 is dominant. To account 
for the Stokes phenomenon we simply replace ^2 with ^ 2 + ^i- Then we change the lower 
limit of the phase-integral to and we replace ^2 with ^2 +^1 again to cross the second 
Stokes line on which ^2 is dominant. Finally we return back to g where 



e l7l V 



ir + (V 712 + e" i712 ) (l + e 2 * 732 ) e~ ir ] *J 



r (*a) 



etna + e -ni2 j e -ir^(fe) 



(5.11) 



The bar on g is to distinguish this solution from ty g . 

In order to return back to e, we reverse our initial steps. We first go back to / by 
crossing the first Stokes line. Then we change the lower limit of the phase-integral to t\ 
and we cross the next Stokes line. Our final solution using the fact that 712 = 732 is 



e -2i 7l2 \ e -iT 



+ i 



ir + (2 + e 2ir 

(l + e 2 * 712 ) e iF + (2 + 2e 2i712 + e" 2 * 712 + e 



Comparing ty e and \l/g requires that 



e 2 * 712 ) e ir 



This will lead to 



2 + 2e 2i712 + e" 2i712 + e 4 * 712 



4*712 ^ g _ir 



(fa) 
2 



(5.12) 



which gives us the correct monodramy in the clockwise direction of the form e * r . 
(I5.13J) gives us the WKB condition 



(5.13) 

(5.14) 
Equation 



,2iT 



-1 - 2 COS 2 7 



12 



(5.15) 



Substituting the values of 712 and Y from Eqs. (|5.4|) and (|5.1U|) respectively, the general 
result is 

(2n + 1)tti + ln[-l - 2 cos((a - 1)tt)] . (5.16) 



Note that this result is independent of the choice of the parameters b and q. However, 
the result is only rigorous for rational b. Nonetheless, this result a common result for 
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Figure 4: A schematic illustration of Stokes (dashed) and anti-Stokes (solid) lines for 
1/2 < b < 1 with q = 2 in the complex y-pl&ne. The hollow circles represent the zeros of Q 2 , 
while the filled circles are the poles at the origin and the event horizon. For u = ujr + iujj, 
with ujr > and ujj < 0, the anti-Stokes line labeled e asymptotes to a straight line at an 
angle arctan(— uji/ujR)/qb and the Stokes line labeled k asymptotes to a straight line at an 
angle [tt/2 + arctan(— u>i/ujR)]/qb. 

the generic class of black holes studied in this paper. As we mentioned earlier, we are 
particularly interested in the physically motivated case where a — 1. For all odd integer a, 
including a = 1, we get the familiar WKB condition 

= (2n + 1)ttz + ln(3) . (5.17) 

J-BH 

For even a, on the other hand, the imaginary part is zero and there appear to be no 
large damping QNMs. It is important to realize in the context of 2-d dilaton gravity, the 
parameter a can in principle take on any non-negative real value, in which case the real 
part of the large damping QNM frequency is not the logarithm of an integer. 

6. A different choice of the parameter q 

So far we have shown how to solve for the highly damped QNM frequencies using the 
change of variable (J4.8|) with the particular choice of q = n, where n is the denominator of 
the rational parameter b. We mentioned earlier that there may be more than one choice 
of the parameter q for any particular value of b, which may produce the desired behavior 
of anti-Stokes lines of the type shown in Fig. El In this section we would like to explore 
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a different choice of the parameter q. If we restrict ourselves to the region 1/2 < b < 1, 
it turns out that with the choice of q = 2 we always have q(l — b) < 1. This means that 
according to Eq. (j4.11|) we always have only one horizon located on the positive real axis 
with no fictitious horizons in the complex y-plane. For this choice of q, according to Eq. 
(I4.12J) . the function Q 2 has four simple zeros. The behavior of the Stokes and anti-Stokes 
lines in this region of b is shown in Fig. UJ Note that the choice of q = 2 is independent of 
the parameter b. In other words, in this range of b we do not need to assume that b is a 
rational number. Therefore, our analysis in the region 1/2 < b < 1 include both rational 
and irrational b. Using the same methodology as in the previous section we can solve for 
the highly damped QNM frequencies, and the result will be the same as in Eq. ()5. 161) . 

7. Summary and discussion 

In Sections 4 and 5 we rigorously showed that for a class of 2-d dilaton gravities the 
highly damped QNM spectrum for single horizon black holes is given by: 

^— = (2n + l)?rz + ln[-l - 2 cos((a - l)vr)] , (7.1) 
Tbh 

where a is the exponent of the coupling h((f)) near the origin, in tortoise coordinates. This 
class of theories includes, among others, the higher dimensional Schwarzschild solutions. 

Although the detailed calculations in Sections 4 and 5 dealt with a specific class of 
models, it is now straightforward to understand the general conditions under which the 
final result ()7.1|) applies. The general form of Q 2 is given in Eq. ()3.4|) . The number of 
anti-Stokes lines emerging from the origin and their angle are determined by the behavior 
of Q 2 near the origin. F — > x 13 , and one is free to choose a coordinate system in which 1—/3 
is an integer greater than 1. It is clear from (|H.4|) that in such a coordinate system there 
will be 2(1 — (3) zeros of Q 2 and hence 2(1-/3) anti-Stokes lines emerging from the origin, 
separated by angles of 27r/2(l — (3). Assuming that F(x) has a simple zero at the horizon 
on the real axis, then the two anti-Stokes lines closest to the real axis can be expected to 
meet on the real axis to the right of the horizon, as in the above figures. As long as there 
are no solutions to F(x) = on the complex plane close to the real axis, the next closest 
anti-Stokes lines will in fact escape to infinity, and the contour that we have used will yield 
precisely the condition in Eq. (|5.1fij) . 

Note that this analysis immediately justifies the conclusions in Refs. [HI] and ^7j, which 
used the monodromy arguments of Motl[10, 11J to calculate the change of phase in tortoise 
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coordinates as the solution circled the origin. These conclusions assume that in tortoise 
coordinates one must use a contour that encircles the origin by an angle of 3ir. This angle, 
however, must be justified in terms of anti-Stokes lines in coordinate x. We now see that 
the relevant anti-Stokes lines in the x-plane are the second ones on either side of the positive 
real axis, and that they are generically separated by an angle of 

AO = 3 , - . (7.2) 
2(1-/5) 1 ; 

When F ~ x 13 the relationship between the tortoise coordinate and our re-coordinate is 
z ~ x 1-13 , so the corresponding angle in tortoise coordinates is (1 — (3)A0 = 3n. Thus the 
results are completely general and valid. 

The general consensus among researchers seems to be that the large damping black hole 
QNM frequency has little to do with the quantum gravitational microstates that account 
for black hole entropy. The results in the present paper seem to support this view: the 
coefficient of the highly damped QNM frequency derives from the exponent that determines 
the x — > zero limit of the matter-gravity coupling, expressed in tortoise coordinates. In 
the general class of 2-d dilaton gravity models this exponent can take on any non-negative 
value, so that the coefficient that relates frequency to temperature is only ln(3) in very 
special cases (i.e. odd integers). For spherically symmetric, higher dimensional black holes, 
the ln(3) seems to be more universal, but again, it is difficult to see how it may be connected 
to the quantum gravitational microstates of the horizon. 

Moreover, if the connection between highly damped QNMs and quantum gravitational 
microstates were real, it should persist in some form for all black holes, not just single 
horizon, asymptotically flat ones. This seems not to be the case. In fact, it seems that 
adding even an infinitesimal amount of charge, for example, so that a second horizon 
appears near the origin, changes qualitatively the nature of the highly damped quasinormal 
modes. An investigation of this mechanism in the general framework presented here is 
currently in progress. 
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